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ABSTRACT 


This thesis develops the theory for and presents a 
digital computer program capablesot determining the natural 
frequencies of a three-dimensiondal piping system having 
Bapitrary conticuratiom., The vamalysismlces stiiesmetnoedmaon 
meamster Matrices. Piping haneers, loops, and conplex 
branches (branches emanating from branches) have not been 
eireluded in the analysis. A distributed mass model 1s used 
for straight pipe sections, while mass is lumped for curved 
mweerions. Inclusion of shear detlection and rotary inertia 
mnopitcional. 

Several piping contreurations are analyzed Usime. the 
program; the resugits ane Sones with analytical solutions 


or values from the literature to demonstrate the accuracy 


peaeintegrity of the program. 
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Loe JNRODUCITG. 


A. THE PROBLEM 

From an engineering standpoint it is important to be 
able to analyze and predict the vibration characteristics 
Sitestructural elements and assemblies. ine mavncnteredal 
aspects of several approaches to such analysis have been 
known for a long time, but the volume of arithmetic opera- 
tions dictated a gross simplification of complicated systems 
Mmonder tO Keduce ther arithmetic to Mandecanlcw miminece 
ane development of the electronic digital computer within 
the last two decades, however, has provided a means for 
tine practical and accurate analysis of quite complicated 
vibrating systems without great simplifications which cast 
abc eon tie inteority of resultSae Minturn. tnisenacmpro. 
vided a stimulus for further developments in theory which 
take advantage of computer capabilities. 

Most methods of analysis of actual systems in which both 
mass and elastic compliance are continuously distributed, 
mvolye a discretization of one Sort or another; usually 
the total mass of the system is considered to be lumped at 
meeinite number of distimet points, tvemiasses being con 
meeted by massless Spring elements. If.is almost universally 
meomed thatsene eclasticsenaractcnis ticceormemesS) Ss tolled ne 
linear and that the excursions of points in the system from 


enerrs CQUILIDrIUM POSITIONS are “SUEPIelentl small. tnae no 





geometric nonlinearitiés are involved... suche anmeidea le zat ie: 
leads to n second order ordinary differential equations. 
Provided vibrations are assumed to oceummi1sochronoucl, 
and no damping is present, these equations reduce to a 
system of n algebraic equations in the meanpleeudecmana 
iPmeeaqucney squared. 

The most convenient and expedient way to cast this 
maoplem 15 in matrix form. Thus “engineerssusine m@nese 
methods deal with such quantities as the mass-matrix, 
PiewStlitness Matrix, the modal matrix eecce. 

One thewother hand, for certain simple cases of dastribuced 
isos and elasticity, it 1S not necessary to lump the mass 
moo a finite number of “eencentrated masses. Instead the 
wo titeron problem Gan be formulated as departialedi teeremmral 
Squatzo0n which can be solved directly. Recently, this point 
Se yiew has been extended for the analysis of more Complicated 
Cases than can be treated by the conventional use of partial 
m@irererential equatvons. Namely, the system beings analyzed 
iedivyided into a number Of sub-systems, a Satisfactory 
amalysis existing for each one; then, by a procedure which 
will be described in greater detail later, these sub-solutions 
are combined to obtain the solution applicable to the original, 
momplicated configuration. This procedure, known as the 
transfer matrix method, evolved from a method described by 
Heetolzer fitty years aco for the@amalysis of torsz1onal 


Vibrations. Applications of Holzer's viewpoint to the analysis 





of flexural systems was made independently by M. A. Prohl 
and N. O. Myklestad about thirty years ago. The essentials 
of these methods were abstracted by a number of engineers 
and mathematicians who developed and formalized them into 
what is now known as the transfer matrix method. The names 
eos, Falk, K. Marguerre, E.G. Pestel faamder. “A. sleckic. 
among others, are associated with this development. A book 
iyaeeesctc! and Leckie {Reft. 9) Prenpresentcsmetcmmnl lected 
most accessible treatment of the subject available today. 

The transfer matrix method is most naturally applicable 
to the analysis of "chain-like" structures composed of 
elements which are adjoined at distinct points and for 
each of which a satisfactory vibration analysis exists. 
ieeesimplest application is to the case where such a charn 
@eeelements has but two ends. However, 1€ can be applied 
©o structures having slightly greater topological complexity, 
DG ewt Simply #senot appropriate fem the analysis of Strae- 
mmmes which are highly reticulated or in which major elements 
are joined along lines or surfaces rather than at distinct 
points. 

Piping systems are substantially chain-like and topologi- 
cally simple. Accordingly, it was natural that the transfer 
matrix method be applied to the analysis of vibration in 
prping. in this country the first such effort which was 
reported was that done by G. E. Fink in his Naval Postgraduate 
School thesis of 1964. At about the same time similar 


efforts were being made in Japan under the sponsorship of 





the Japan Society of Mechanical Engineérs. ft 15) known 
that the Japanese have produced useful digital computer 
programs for the analysis of piping using the transfer 
matrix method, but it has not been possible to learn any 
of the details since these programs remain under the pro- 
peletanry control of the Japanese cevernmenc. 

Fink developed a program whieh has "eeonveniemerccn 1a) bey 
mma olves good results, but which 1S limited to systems 
iving entirely if a single plame. His prograieis capable oF 
making separate determinations of "in-plane" and "out-of- 
plane’ vibrations, which, for the case of uniplanar cen- 
figurations, are uncoupled. Further werk on the theory 
Poromecone by W, S, Barrd, Jmeandyd) Ly simmons im) thers 
Naval Postgraduate School theses; this work mainly related 
messystems having appreciable topological complexity.” In 
ims Naval Postgraduate School thesis, Y. S. Kim returned 
to the development of Fink's program, adding several useful 
meatures, including the incorporation of several alternative 
iaenematitcal techniques for accurate solution... Some ot ‘tliese 
techniques reduced solution time below that required in Fink's 
oracginal program; however, Kim's program is still limited to 
uniplanar systems. 

Mecording Ivete Writer Under COO kmromce velo mca Oo am 
enable Of tréeatine truly three-dimensienal cConticurations. 
The remainder of this thesis presents the analytical back- 


sround of this development, lists the program which includes 





rules for itS operation, and gives evidence sconcenninc strc 
Nileh ~accuracy .O1 she vesulece 

It should be made clear that no claims are made concerning 
time relative usefulness Gf EheG transter macrieemioenedmas 
esompared to llumped-mass methods whweleamespuese mit) mann 
widespread use for analysis of piping vibrations. In 
memeral, the writen believes that the nestults when eare 
Optained by the transfer matrix procedure are substantpally 
more accurate than can be obtained in most cases using lumped- 
mass analysis. However, the delicacy of the calculations 
mivoived in the transfer matrix method seems to inply 
mimait Only a limited number of the lewest frequencies (and 
Met: associated modes) can be Obtained usine present pro 
gramming methods and computer hardware. Greater word-length 
meeene computer, or availability of fast multiple precision 
capability, would permit obtaining a larger number of 
meecquencies and modes via the transfer matrix method: 
mover iy, the motivation for the development reported 
mocim has been simply that the transter matrix method 1s 
elearly appropriate for the analysis of piping vibrations 
piemchatethe application should@he made. it seems biakely 
miat in some cases and for some purposes the transter matrix 
method may be more convenient and economical while in other 


eeses lumped-mass methods may be prererabic, 


0 





B. THE TRANSFER MATRIX METHOD 

opecific details of the applications. trans -cen Maen. 
methodology to the analysis of piping vibrations will be 
given in subsequent parts of this thesis. Background 
material and a general exposition of the method is most 
readily available in the book by Pestel and Leckie [Ref. 9] 
meade it would be pointless to ineélude Nerecinead parapmnase 
of such material. However, for the reader who is not familiar 
teen Che ideas behind the transfer Matviseemethod, ay veny 
brief exposition is given in this section. 

MileapprOopriate “StauLe vector aS dered winvtenmo1esome- 
force-type and deflection-type information capable of specify- 
itmche state" at each of several points in the configura- 
tion at which sub-elements are joined together. These 
imamtities Specify the configuration and internal force 
system in the structure at its extreme isochronous, deflected 
Mesition, For example, in a simple torsional system, there 
m@mer,Only two quantities in the state vector, the angle cf 
moration £rom the equilibrium position and the terque 
MEiismacted thom one Clement wre the next... ihe elastic 
fomglrance and the inertial characteristics of an element 
Bermit obtaining equations relating the quantities in the 
meate vector at the left end of the element to those in 
the state vector at the right end of the element. When 
these relations are put in matrix form it is found that 
the state vector at the right can be written as the product 


Cimene State vector at the lekt premuleinitedsn, aesquane 


ES 





matrix, called a transfer matrix, the elements of which 
represent the elastic and inertial properties of the element. 

In applying the method, a state vector is chosen for the 
left end of the assembly of elements in such a way as to 
satisfy conditions of restraint at the left. Then, by suc- 
cessive multiplication by transfer matrices representing 
individual elements encountered in order, progressing from 
the left end of the assembly toward the right, one arrives 
mame representation for the stave Vector ede then meneiemcid. 
The individual transfer matrices are functions of isochronous 
MmenMcncy: accordingly the rrehe endematryx seas sinet nonmon 
maeauency. When frequency iS adjusted, any particular 
meequency which permits satisfaction of restraint at the 
right end 1s a natural frequency of the system. The several 
meme VECLOrS at the junction points ceive the deflected 
Soerouration, or mode shape, corresponding to this trequene.s 
when this value is substituted in the expressions for the 
meaue Vectors. 

Ome Simple cCaSes it 1S) possible to weary somt cient letee 
Wimeattons in literal form, keeping the trequency jas an 
unknown parameter. Then the right end restraint conditions 
provide a polynomial equation in the frequency which can 
be solved by customary methods. However, in relatively 
Complicated cases it 1S not feasible to retain the frequency 
mised literal parameter. in these easesua derinite mumerical 
value is assumed for frequency and this leads to a numerical 


measure of the failure to Sadtasty pucnteenas res tlaint 


lee 





conditions. By assuming different values of the frequency 
and by constructing a curve of measure-of-failure versus 
frequency, the correct frequencies seanm ee mop -oumcae 

The transfer matrices which describe the individwal 
elements may be a theoretically complete and full repre- 
sentation, or they may themselves by approximate representa- 
tions. In the application to piping systems it is convenient 
Mmemrepresent a straight léneth of pipewby sa teancecr marie 
rarely 1S theoretically complete’ and exacty § For Ssysecic 
Somposed of only straight lengths of pipe there 1s ne 
approximation whatsoever in the transfer matrix method 
former than that implied by internal round-off in the com- 
puter); this may be compared with the fact that all lumped- 
mass procedures definitely imply an approximation and thus 
mesult in some theoretical error. 

Curved pipe elements (elbows and bends) could be treated 
meeaGily Since adequate theory exists for the: construction 
eo: theoretically exact transfer matrices with which to 
maoroseit  chicm. However, | tniemwema be = Se vetCmConpl bau Lom am 
ier reculties involved with such exact representation. 
Accordingly, with the awareness that the total length of 
eurved pipe 1s but a small fraction of the total length 
Se both straight and curved pipe in the average configura- 
Mone lt was decided to Pepresent Curved pipe oy aussucces— 
Sion of massless curved sections having proper elastic 
properties, alternating with appropriately located point 


masses. For each of these idealized elements, it is 
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possible to obtain the corresponding trans temeniaentec om oar 
no computational difficulties. Thus, elbows and bends are 
represented in much the same way that as edome wwii tiesenc 
iumped-mass methods mentioned prevveus |e ets be luciyca 
that the approximation introduced by this point-mass 
mepresentation of a small part o£ a Contigurationsisequite 
moncraple for practical engimecrane punpescs -—meloweveT m1 
mecettainly possible to modighy the programe neportcasmenein 
memas tO provide for exact representation Or GUnyed selc — 
Ments. This would not involve major changes but only the 
fmevelopment of a subroutine to replace those given herein 


for curved elements. 
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IIT. PROGRAM DEVELOPMENT 


A. THE STATE VECTOR, TRANSFER MATRIX, AND FREQUENCY 
DETERMINANT 


Pe ouacte VeGron 
ies Sitamcmuce £O 1 z2. at pOInt Deer aneelast iets > com 
ma column vector whose components are the generalized 
frsplacements , D., and the generalized forces, F attra 


Mount, iin matrix notation, we have 


For a beam free to move in three dimensions with 
Monmeltudinal, torsz1onal and flexural elasticity, the state 
vector takes the form 


- deflection in x-direction 
- deflection in y-direction 
- deflection in z-direction 
ge es\Balelgm ehoveyje X< ~ eo cLS 
rotation about y-axis 

= roOtattOneabouc z-axis 


=" 9X14] Shores Inex-ai1receron 


<_ 42 oee = «= CF 
i 


= Shean LOree=iiey -alrcetLon 


~< 


="Shearecouvec iimez-cilnee tron 


N 


Torque about x-axis 


- Moment about y-axis 


cx lel 
— 
1 


Moment about z-axis 


= 
=. 
N 
I 


Mme erder in which the individual displacements and» 1o1rees 
meeur in the state Vector 1S /0r NO parercul an wmporeanee 


so long aS consistency 1S maintained. Program VIBREL 


les 





(for VIBRation of ELastic piping), the program developed 
herein, employs the arrangement given above. 
Z. Sign Convention andetcordinatess aaem 
A right-handed cartesian coordinate system is used 
to describe the local generalized forces and displacements, 
Maemx-axis coinciding with the centroidaly axis voters cmpa pe 


gemeontown in Fig. 2.1. 





Vz 


bigure 2.1: Vectorial Representation ot State Vector 
Displacement andwrerwees 


AYcut Made ap any locaton sot che pape wil lex posew eo 
mae@es; the face with its outward normalvinm the positive x-axis 
m@imecetiOn 1S Considered positive. The state vector parameters 
EygempoSitive if their vector representations coincide with 
miceairections prescribed by the cartesian coordinate 
fyotem at the positive face as) in Fig. 2-1]. Rotation and 
moment vectors are portrayed according to the right-hand- 
screw rule by a double line with an arrowhead. 

Se lranstemeiatrvecs 

The transfer Matrix. 1s ede imed ea scepeincmia tod evneecea 


melates two state Vectors @t SUCCeSsiVvenselleetca poImresmon 


16° 


division in the elastic system. 


ieee Dera 
it 1 hall 


where U. is the transfer matrix... Uhewtyanstcr  ilacniecm ct 
be square, and in the case of a spatial system must have 12 
rows and 12 columns. When the transfer matrix relates the 
mlexibility between the points a7) and 9-1 cf auconeumugu. 
Poem it Will hereafter be referred tovas autie dicen wen 
mnen a discontinuity arises in either force or displacement 
as with a point mass in a lumped system, the transfer matrix 
felatineg the state vectors on either side of the disconmting- 
[igents Called a point matrix, Po) Ine term “transtenr smiaunaen 
vee henceforth be used with regard to ™em@e matrix, Uy whieh 
imeludes mass and flexibility in the rellationship between 
fre “state vectors at two distinct points, A complete discus. 
Sion of transfer matrices, with examples and derivations, can 
Memround in Ref. 9. 

4. The Elimination Process and Frequency Determinant 

The process of combining the system transfer matrices 

moetind natural frequencies can be demonstrated by a Simple 
planar example using the Myklestad method of replacing a 
distributed mass beam with lumped masses connected by flexi- 


bee elements (Fig. 2.2). Assuming that the field and point 


0 1 2 
Lor 
ES ™ ET 
Figure 2.2: Lumped Mass Rennes cimeatyommou 


a Planar esean 


ee 





matrices for a massless flexible beamsanawaspemmemnacoma ne 


already known, then 


L Re ee eee 
24 Viz» Zz, = PZ4; and Z5 = V524 
W 
where z2. = Y and. superscripts b and Ro andseatemlcteucide 
ie 
M 
bs ge 


Hierioghnt Side, respectively. Combining thevaboveryields tite 


hmelationship 


ee a 


Z» oPVi2, = Uz, (2-1) 


tere i 1s the system transfer matrix. Written an full 


Dies relationship 1s 


il Sl, > 126 71 5eaenel " 
Y aie 22) 72 ae 
al) | ee Maier tars Sg | | Ys 
: Jalen aie ei) || 


sa | 92 
the ae being vibration functions of mass, stiffness, and 
meequency. Application of the boundary conditions B= Wy = 


w 


a > = reduces  themlatrixeedlat lone com tie ssell own 


= (2-2) 


Because the elements of this square matrix are functions of 
mmeawency, its determinant 1S called tie tre@mency detemm- 


nant of the system. Note that this 2 X 2 matrix is a 


18 





submatrix of the 4 X 4 matrix appearing above. A non- 
trivial solution of the equations 2-2 regimes Ehaeethe 2 
of the frequency determinant (A) equal zenro A vale cree 
frequency, w, which causes A to vanish is called a natural 
PmECGUCHC Vy wOL, Lies s tiene 

The natural frequencies could beusemicd eran edimece mi, 
mEovaded the expressions) for, tae 255 Werte, Kept in) temo ou 
mme circular frequency, w, but ingpractie@egenesc expression: 
aoe forms much too unwieldy to bewiamawed) ciate ti). 
Program VIBREL evaluates the frequency determinant of 
order 6 X 6 for successive numerical values of frequency, 
Momo values of the frequency determinant Corresponding te 
maiculral frequencies of “the system. 

letethe state vector be Composed ome Beauamt) Gco me eae 
memeompute U ain Eq. 2-1, two multieiications of Zr Moan 
ieyerices aré required. Noticing that application of the 
boundary conditions to the state veetor at point Yewi el 
eause r of its elements to be zero, we find that r columns 
of U are multiplied by zero, and thus play no part in the 
werculation. Kim [Ref se6))| imtrodiuced  thesconcep a on gene 
"state matrix," S, such that the left end state vector, 2, 
meer quantities can be written in terms of a compressed 


Beate Vector, ae Nav ie wonder ecuamit ist ics amy ez 


meme S 1S thesstate matrix at thie, lenemend. Nayving 2 


columns and 2r rows, iltSmelementsmconstoeingmenmis got Zeho> 


i 





and ones. All elements of the state matrix are zero except 
that in each column, unity appears in the same row position 
as the corresponding non-zero term in the state vector, 
boundary conditions having been applied. If the state 
matrix is used in the multiplication scheme, then we can 
modtiply zr X r matrices by 2r X Zr matrrvees and reduce 

the calculation considerably. The resulting 2r X r system 
Enanstfer matrix requires only the applicagaon of right end 
boundary conditions to yield the frequency condition. 

Teo illustrate this state matrix CcConeepmemest Sinply, 
consider again a two-dimensional case in which the state 
fmoetor consists of four quantities. Them three dimensional 
case is an obvious generalization of the following discus- 


Sion. Taking strictly the state vector approach we can 


j 





Figure 2.7555 - si nelle sha ane heer deaaidien Sy Siem 


express the relation between the state vector at Zz and 


Mmirac at 0 by 


25 = U,P,UZ, (2-3) 
where U, and Us are straight section transfer matrices and 
the point matrix Ps accounts for the effects of the branch. 


Carrying out the matrix multiplication and exhibiting the 
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System transfer matrix and state Veetorsmywiichmeconsacden 


bending only, we have 


" Slit ile Sos “Se 0 
v 254 422 223 474 0 
= (2-4) 
UE 437 *32 455 234 Ay 
M , [4a 242 9843 344 Ao 


The boundary conditions have been applied to the left end 
meee Vector; the quantities Ay and A, Gepreseme, them emcee 
Gmantities of VA and M, Trespeceively —MeApp lea Glo Me eesti 
right end boundary conditions yields a frequency condition 
of the same form as Eq. 2-2. 

Heime, the state matrix approach anctecaqemrd a7. 
becomes 


= * = 
Ae WP Se. as 


where S is the state matrix z” contains the non- 


0 
0 
il 
0 


-= Oo O&O & 


moeeergquantities or the left end State vector. AREer multi 
Mmmeation of thesmatrices of Eq. 2-5,8the number of andivaduae 
mieltaiplications being considerably léss than before, the 


mollowing relation is obtained: 


a Ss 14 Ay 
v 3 ns A 
= |, : (2-6) 
a 
Wz 33 “34 





We see from this result that inclusion of the state 
Matrix compresses the Left end State Vector toe. duane 
while providing the same useful information as the matrix 
on Eq. 2-4, inasmuch as €he same frequcncey condi wens 
obtained after the right end boundary conditions are 
applied, Also, if the relative magnitudes of the unknowns 
A, and A, have been determined, the mode deflections and 
forces at a point can be calculated simply by multiplica- 
muon Ot: the compressed left end state vector, consisting 
Pie diantities, by the Zr X r SsyStem transter matrix up 
momthat point. The appropriate mode frequency 1s used in 
mie computation. This will] be discussed in greater detail 


miechapter 3. 


fee -) PLPING GEOMETRY 
1. Working Points and Piping Nomenclature 

The geometry of the piping system as it is used in 
program VIBREL is based on working points. A working point 
iemnene System occurs where cross section or properties 
change, where branches join the main member, where the 
Buping changes direction, and at the ends of the main member 
pma@ Dran¢ches. These points are measured in cartesian coordi- 
ieecs cron an Orlpin Pocated@arbitrarity but fixed in space. 
For purposes of this discussion we shall distinguish between 
mre two types Of working points at which piping direction 
changes. At a corner point (Fig. 2-4) the piping direction 
Ehanges with no discernible radius Of Curvature. At a 


Deamon htO sees), Which 1S denineddsmunemimtersect rom 


Diem 





of the tangents to the two ends of a curved pipe, the radius 


Of —CUTVature 1s = inibece 





CORNER BENO 
‘ POINT POINT 
meoure 2.4: Piping Corner FP OUMe 2 peo INO ecm 


The main member is designated as the longest continuous 
run of pipe from which all branches emanate. Program 
VIBREL cannot handle branches emanating from branches. 

Z. Unit Vectors and Coordinate Transformations 

iva set on mutually or thoconalex,.) wands 2a tmd e 
meetOrs are formed at the main member and branch Starting 
mepeesS, and prior to and following every every directional 
Giange, we can determine, by angular differences in succes- 
merece setS Of UN1t vectors, the appropriate coordinate trans- 
imomiiation angles. the local coordinate systems described 
by these unit vectors are formed with the piping system 
miei1ts quiescent configuration. 

Let us take, for example, the two sections of pipe 
pictured in Fig. 2.6 determined by the three working points 
tweet, Gi,where F 1S a corner point. Also let the two sections 


of pipe, EF anda FG, be represented by vectors K and B, the 


(6) 








Beeure 2.0: Vector Representacivon Omelane et mom aaa ec 


Samponents Of whichmehaye been GeLermined a1 Onesune coor ai 
Mees Of FE, F, wand Gow then theswanele 7 can be derermmimedm ry, 


the following relation: 


bo 
1 


ly] = Cos” 


|B] [Al 

The sign will be determined later. Since the vectorial 
representation Orectadte Vector eGU dir iee die cmmin ns 160 remo 
maeuires that the x-axis coincide with the centroidal 


fears Of the pipe, let 


Hem the purposes of Succedin@ derivations we desire that 


ey, be in the plane of A and B and point outward from the 


g@irectional change. Hence 


oe aX =. wows 
et ae and e. = ee alee 
EL ea) Y 
By the right hand screw convention the rotation of amount 


fe will always be in the negative ce, Senses. lhe: set Of sUunit 
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VeClLoilae Oe Uy» U, are found simply by rotating the initial 


unit vectors about the locaimez-axis by —ehemane leu 


Now consider the ditfterence in orientations betveecn 


the Us U,, » U, vectors and those prior to the next directional 


change farther down the pipe, U,,, Uys s and Us. These have been 


constructed in a manner similar to ey ey» and e.- Referring 


to Fig. 2.7, the angular difference between the 





Figure 2.7: Unit Vector Rotation about X-axis 


moomadjacent Sets of unit vectors 1s merely the rotation @ 
about the local x-axis. Alpha may be determined by the 
expression 

ee 


= Cos. ae 
| 0. | os ~ (UU, 


”) 
iis rotation 1S in the positive (negative) om sense accord- 
ing to whether the triple product (U, x0, 1) UL! is positive 
(negative). 

In this way, the state vector is always expressed 
miethe appropriate local coordinate system. For the simple 


eyscem Of Fig. 2.8 


LS te 








Figume 2,6: Planar System with aneeurcmalmoeanac 


iN 
eS 
N 
N 

be 

i 
(=n 
N 


which combining gives 


Z5 = UT 


where UL TSeecrew trans tOrMat 1OlImiarn Ix ec Orrespendanm: 


rotation y about the z-axis. For a 3-D system with 


moeror detined aS in Section Z2-l this matrix 1s 


Le = DIAG Eee Nios wy vn ; where M 


Emmibarly, for the x-axis rotations 


tae DIAG [T_, T., To, 1} , where 7, 


ee -D TRANSFER MATRICES 


COSiy soln 
Syn ay 


it 
0 
0 


ime otics senate Pipe sivans ier iia tama 


0 


0 
COS a 
-SIN a 


COS 


0 


eOmd 


the state 


In order to maximize accuracy, a distributed mass 


model was decided upon for straight pieces of pipe. 


The 


peocram wses the straight bar transtem matrix weited in tehne 


Poeddi sOiecr. 9. bDUt revised £6 conform to the ordering 
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of the State vector described herein. These) ised aise 

which appears in Appendix B, combines the vibration func- 

tions for longitudinal, torsional, and flexural modes and 

Gan. include or neglect rotary inertiasand ss meat sden lee auome 
Z.. Develeument of the Curved Pipe. bme Vda tiie. 

Due to the complexity and increased machine time 
involved in consideration of a distributed mass model for 
curved sections of pipe, a lumped mass treatment was chosen. 
This, coupled with the distributed mass approach for straight 
pipe, yields staisfactory accuracy combined with minimal 
Somputer time. 

Given the compliance of a piece of massless curved 
pape at its center of curvature, we seek the field matrix 


te ae 


Weouch that z, = Vz, (see Fig. 2.9) where z, = { 


PA 
4 N 
/ \ 
/ \ 
/ \ 
a ‘ 
/ \ 
/ \ oy 
y, " 


Figure 2.9: Maseabess @@unved™= Pine 


eompliance matrix a is defined by the relationship 

a = co where uh iS a vector of generalized displacements. 
observed at P and ES 1s some vector of generalized forces - 
ool tcdeat web ountc P 1S assumed POmecmietdencauhy wc reat 


Pe 





bar to one end of thespipeswnilte the other vende encom pap 
remains fixed, The compltance matrix Lord) pipe moamanua 
mes derivation.are presented in 2 paper) bye) ee eeemock 
Ree. 2]. 

The analysis requires that we findweneweonslranee 
maerix at point L which 1s obtained by the appropmiaece 


fgeteatioOn and translation of coordinates from point P such 


fina t 
- T 7 
Cy = Boy Lpy Cp Lpy Bpy 
I 0 0 Zee Viren, 
where B eG = Poul LP 
PL Dor I PL Ze = 0 ya 
Dae Pe POL 


COS (xp x7) COS(xp,y,) COS (xp ,2,) 
Loy = » Kp, = |COS(yp,x,) COS(yp,y,) COS(yp,2,) 
COS (z,,x,) COS (Zp y,) COS (zp,2Z,) 


ite Matrix Cquation for the seneralized displacements at 


meint L can now be written. 


: T 
DD = CFL * Ber Lip Pp 


winere the rotation and translation Matrices use the above 


MOotcats1on, Solving Lor Dy we arrive at 


Com 





a eee 
De = Lip Bre Py - bir Burr Sb, Fr (2-7) 


the Cquation for the generalized forces ae spice woe 


Peete. pape hasuet memeno ten 


FR = ~Lor Bor oT ‘ CS ss) 


If equations 2-7 and 2-8 are combined in matrix form, we 


have the result 


T | T, T 
Dp Sone ier Mik 
| eee ba ------ 
= 
He t enn oT 


Pyaeteie Sion ConventloOnedmseussed in Section: 2.) point - 
is located at a negative face, so after appropriate sign 
Gmamoces are taken into account, the field matrix for the 


bend becomes 


lee aa Taal 

op cas 0 Ss Bas SO 

Vs |e-------b--------- 
0 | ~“RLPRL 


ims matrix Jn its explicit form, as used in program 
hebREL, appears in Appendix B. 

See weve LOplenieOouecunved= line POnuntetia tr him 

The mass of a pipe bend or elbow is taken into 

meeount through the transter Matrix pmueeedure by Jumping it 
at the center of mass and including its effect on the system 
PeapotiewO (Flow 2.10), the vleftrcndwom elem elDoWw ese ecrlon. 
men resulecing poine Matrix relates Jonlyaune torees across 


pormte.O since the deflections are Continous, 


Xs) 








eee ek) eer Curved spate Showing Location 


of Center of Mass 
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A force sSumilatroneae point 0) yields thes tomiowais 
scalar equations (the quantities appearing in this develop- 
ment are defined in Appendix A, subscripts L and R refer 


me the left and Tight of Point 0, respective ae 


2 : 

Np = Ny - mw (u - y8) (2-9) 
Vi = 

Vor = VoL -mw (v + x@) (2-10) 
2 ; : 

VOR = Yor -mw (w - yd - xw) C2 = 11) 


Similarly, a moment summation at the same point produces 


z 7 2 : 
oy = I, wo (2-12) 
Mop = My » 1 ww (2-13) 
ees Ce BR: es 
ZR ZAG ZO 
Compining equations 2-9 through 2-]4 
fy = 
where 
= -muw 0 0 0 0 -mw “y 
0 Tne 0 0 0 aoa 
0 0 ivi muy mu X 0 
0 0 0 ae eae 0 0 
XO 
0 0 0 0 ie 0 
VEO 
2 
0 0 0 0 0 =a | 


oul 





and. the resulting point onan cauttregeic. 


An inertia tensor analysis was used in the determination 
Seeche Mass moments Of Sinertia about pomnus Come eeean 


meacily be shown that the imertia tensor at point © as 


—~— 
— 


-2 chy eee 
Mee See NY ae 271 e | Tae UNG ye 


BOAR [y-SINy]u + pAR”[ (-% + SIN ie 


Z 
Mea-cosy - 22%) (ij+ji) - Yt - asiny + 247433) 


where it is assumed that pipe wall thickness is small 


sempared to pipe radius. (fhe symbols used here are defined 


meenppendix A.) Then Lees ale 


T 
mdr = k « | © ki where i, 7 
ZO O 


ia La ae 
* Il, Lio = J I, Js 


, and k are unit vectors 


O 


peeome) the local coordinate axes at point OO.) The ensuine 
Mements of inertia and point matrix are presented in their 


entirety in Appendix B. 


D. BRANCHED SYSTEMS 

Any branch joining a piping system has a significant 
emeeecct On its dynamic behavior. While the detleeétions 
across the junction point are continuous, there is a dis- 
continuity in the forces, the magnitude of which is 
dependent on the displacement at the junction point and 


~me nature of the branch. 


2 





Consider the simple biranehe system of eee 





\ Se \/ Uz" 
Se B U Y 
A =a 1G 
a. | 
Uz 
U. 
P U, 
E 
D 
regore 2.11: Unit Vector Orientation at a Branch Junction 


Using the state matrix concept discussed in section 2.1 we 
have the following relationship between the state vector at 
point D and the state vector at point B (the bar denotes 
branch coordinate scheme). 

12 X 1 6 X 6 


B it = Se GOs 
=| -e-x [lag i 
B , 


D 


E 


The zero quantities of the state vector at D have been 
eliminated by application of the boundary conditions and 
imeroduction of the state matrix, and the branch transfer 
matrix has been partitioned into two submatrices Ry and Po. 
The local coordinate system of the branch at point B is 
mepresented by the unprimed set of unit vectors. From the 
relationship 2-15 it can be seen that D. = Ry 2p and 


Fp = Ro Zn» 50 that when Z1y is eliminated, the equation 


between the forces and displacements at point B in the 


unprimed coordinate system becomes F, = Peas (2-16). 
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The vectors of thevgeneralized. forces and displ acenenme 
in the branch system must now be transferred to the main 
member svstem represented by the unit vectors s in Fig. 2.11. 
This can be achieved by one z-axis and two x-axis transtor- 
mations of the unprimed coordinates. (Ne wtivst seonmercnc 
Boamed set of unlit vectors in a manner Identical comenaee cor 
BGarner point (section 2.2), the three working poimes 
feceribine the corner being E, B and ©. gan x-axis) rotation 
oT will align the unprimed with the primed coordinates. 

Next, a rotation y about the z -axis will allow Uys to 
eermmcide with the centoidal axis of section BC resulting in 
mie couble-primed coordinate scheme. Finally, rotation 
mmnouch the angle a about the x'-axis will align the branch 
System with the main member coordinates, s. The total trans- 


mormation Of the displacements can be expressed by 


Dp = G, G Su, Dp = GD, (2-17) 
1 0 0 il 0 0 
where Gy = DIAG UM COSd Sao tNa |r Se0Sa wacinc : 
OR SG ee Ooo east Naw Gose 
J 
eamilarly for Say? and 
COsc SINoeae COSsapesl Na 0 
e SS DEAGHIooUNG €COSco Cs ieee Osc 50 
0 0 1 0 0 Il 


L 


mexewise for the forces, a = GF, (on cary ano ine 


relations 2-17 and 2-18 and substituting in 2-16 yields 
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= Il _ =) = ll 
G Fe = RoR, G “Dp 
Noting that et = gt we now form the force displacement 
equation 
. sak AL ail 
Fp - GRR Gl Gyeaee 


Referring to the free body diaorameain hie eZee ees! 


S B 
a ee 
<— Pek 


fy > F, 


Bemnre 2.12: Generalized Forces at a Branch Junction 


mie the main member force equation at poimt B can ble written 


R L 


ee ee 
Lee py 
= Fo + GRJR G41 eye? B 
itis. Since Dp = De = Dae the point matrix relating the State 


meerors to the left and right of a branch junction point of 
a 3-D system 1S given by 
= fanaa] 
T* [eeaaitetyotet) 7] 
The foregoing development assumes that there is no 
emvature of the main member or Dranch at point B and that 


mo abrupt change of direction of the marn member occurs at 


that immediate location. 


os 





For two or more branches joining the main member at 
a single point, the system transter matrix Us eo r porno 
the branch point is Simply multiplied in turn by the point 


meberix for Gach of the branches. 


pee COMPUTER IMPLEMENTATION OF TRANSPERSMARRIGES 

It would be anvextremely tedvoussyjop pro attemperc.p ! mea 
mermolace ot the elements of the systemetriancicremarnrix wim 
moms Of the circular frequency. SubSeaqwent solution for the 
meets Of the expanded frequency determinant would also be 
meractical for all but the simplest of systems. It is 
men these reasons that digital computation becomés a 
necessity. 

Puoctam VUBREL. which Incerpouates = thes tnceny son 
fieeccaing Sections, was developed to utilize the speed with 
fmten the complete system transfer matrix can be formed time 
fewer time for various values of frequency. A typical graph 
Peeercaucncy determinant A(w) versus circular frequency w& 


meponown in Fig, 2.13. The values of w at which the graph 


A (w) 







NATURAL 


y FREQUENCIES 


fe 





\ 





Figure 2.13: Graph of Frequency Determinant vs. Frequency 
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crosses the horizontal axis)» mamely those values vor stom 
which A(w) = 0, are the natural frequencies of the system. 

The program, presented in Appendices CG ythrouch er sae 
eeded using Fortran gdancuage (Ret: 7 ]\ inedoublewmcenc mam 
moeminimize the effects of round-ecff and) tuumest onecnmeuce 
pamce the prime objective was to produce agworking program 
Semcaemonstrable accuracy and reliability, limited attention 
was focused on the niceties of programming intended to 
meamce execution time. However, Kim's state matrix concept 
mermcluded and this effects a substantial saving of time 
@secompared tomjother procedures. There 1S no reason to 
femneve that the actual programming is particularly ineffi- 
mmemt, however, it 1S not unlikely that close serutiny of 
procnmam details might point out some places where slight 
mererolons would effect economies. 

ihe transfer matrices used in computation of system 
meeauencies are those presented in Appendix B. 

input formats gwere deSigned to provide safistacrory 
weet rcity for the user in preparing data decks, while 
Omenit formats were chosen to check for correct input as well 
as to provide a detailed summary of results. The I/0 


mommats are discussed more thoroughly in Appendix E. 


Sy) 





Cit DLS CUSS LON 


hee SOGOPE OF THE Sse Lie 

In its present form, program V1 BREI=rcecanabn lomo. 
analyzing a piping system having few restrictions regarding 
topology. The system must, however, be composed of a 
main member from which can emanate as many as fifty branches. 
The program as listed in Appendix E can accommodate no more 
than two branches joining the main member at a particular 
pemmt, In addition, there can be no branches emanating 
meen Other branches, nor can there by any curvature of the 
fan member or branch at a branch junction point. 

The total number of sections capable of being analyzed, 
Meesently set at one hundred, can be increased by foilow- 
ieee che procedure described in Appendix €. Likewise the 
member Of branches can be augmented from its present limit 
mamtity. For each increase of one section, eighty addi- 
tional bytes of computer storage are needed above the 140,000 
igeee> New required for the entire program. For most computer 
installations this allows considerable leeway in the size 
of the piping system which can be handled by the program. 

Louprovidesan tdea of the computer tame involved in 
fence frequency analysis, the twelve-section piping system 
Sumcie Sample problem in Appendix G required about four 
minutes for computation of four modes and almost fifteen 
minutes for eighteen modes. The IBM 360 computer compiles 


Pies program in thirty-nine seconds. 
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The usual assumptions of linearity as well as homogeneity 
and isotropy of the structural material have pecnerieemponaecd 
in the analysis. All sections of the system are eassumecamie 


vibrate isochronously with negligible damping. 


B. NATURE OF ERRORS AND INACCURACIES 
Ie Jounal hase Isiciee aes 

AS piping System complexity In@weacc cme sui 
eammultiplications required to form the system transfer 
Meaecix grows proportionately. Each time a multiplication 
occurs, there is some round-off and loss of accuracy; this 
[matunction of the Sienitiednt figure Capacity Of Ghe comouver- 

FORTRAN, in conjunction with the IBM 360 and double 
Mrecision arithmetic, has the capability of carrying numbers 


Suelo Significant digits with an é€xponent range of Troma to 


ne Simpler SVSvenmsvecane pe Checked Poy sAaGeutdey saincmrinimaic 
Eemeacnt that this significant digit capacity keeps round- 
Gmemerrors negligible. Since with more complex systems 
eeler forms of errors creep into the results, we have no 
foeeeot predicting the exact contribution of round-off 
emeors tO total inaccuracy of frequency. 

It is reasonable to assume that round-off plays 
aemore significant role in the 3-D piping system accuracy than 
in a planar case with equivalent sections because of the 


(meater number of coordinate transformations that are neces- 


Banye tOr System description. 
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2. Lumped Mass Idealization 
For sections of curved pipe, lumped mass connected 
by massless springs 18 used as a model hor the reales, seen 
Haas idealization leads to inherent inaceur acy sinernomcmae. 
Sis of systems containing bends. 
oe 2eros 0% the Preauency Vetermuname 
Oneée it has been established that the frequency, 
determinant has crossed the zero axis, Program VIBREL uses 
Ma@etier's method of successive bisection and inverse parabolic 
immeecr polation £0 détermine the zero. Care should be exer- 
eed In considering as significant only the same number of 
feeaees aS in the acceptability criterion which was specified 
mmmecne SOlution. (See Appendix C.) 
weeovcncr soukees Of Imaccunpacies 
THe wImperfect nature of Lhe =reaiesys penmeandsimie) eu apm 
eamens in the measurement of working points will render the 
computer solution only a good engineering approximation of 
time can be expected in actuality. In the discussion of 
deeuracy 1m Succeeding sections, inaccuracies caused by real 


Sevem imperfections will not be considered. 


SeeeeaCCURACY AND INTEGRITY OF SOLUTION 
Pee SOLUE VON AC CUTraCcY 
COmetdenGesine thes ACelmiel mG ReenomOnOer ales OllGilon 
Pereestablished by three comparison tests, the results of 


mien were tabulated in Appendix F. 
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a. Analytical Comparison 

For single straight sections #€ pipe a closed 
form solution was obtained from the governing fourth order 
differential equation with the appropriate fixed, free, or 
ball joint end conditions. The largest variation between 
bie natural frequencies calculated analytmeally and those 
computed by program VIBREL was 0.009%. 

Bb. Retememecer Value Conga som 

Where the exact closed form solution was not 
availiable, recourse to the literatume prowided some compari- 
son frequencies. For the simple curved pipe shown in 
Memendix F, Refs. 5, 6, and 8 contained a range of frequency 
moameacs for the first two medes. VIBREL «results for the 
game curved pipe fell within this range in each of the two 
modes examined, with an average variation from the reference 
Peeees Of about 4%. Because, for a curved pipe, program 
feoreL includes bend flexibility factor and an accurate 
memene £or shear distribution factor taken from Ref Z, the 
values computed are considered to be at least as accurate as 
fmese contained in the literature. 

For a SingWVembranch System composedmot stmameht 
pipes only, comparison frequencies for three modes were 
faeeelable from Ref. 6. The largest difference between 
Pererence and VIBREL in this case was 0.06%. 

Cc. Dual@inalysiseConparason 
Dual analysis is based on the fact that the system 


Eeamster matrix is directionally dependents One end of the 
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system main member is designated as the starting point; then 
the component transter matrices are multiplied tecernier san 
Sequence proceeding toward the other end. It angennon 
were to occur in the structure of one of ‘the transfer 
Mmaeraces or an the méthod of their combitiatvonmeomeon er 
eworem transfer matrix,the results frompstartineuaeeorppes ine 
emas could differ significantly. 
An improperly constructed state matrix would 
also be evident from a dual analysis. 
the straight, pipessys tems andes namemeod sis tems 
ehowed no Significant differences when subjected to dual 
analysis. The maximum difference observed in eighteen 
Medes Of the complex system dual analysis was 2.6732, 
while the average difference was 0.58%. 
ee oOGluctOn Meer rey 

Solution integrity was established by observing that 
the mode frequencies detected bv the program were, in fact, 
those of the system with no omissions. This was accomplished 
eyaecomparison with analytical and reference results. Dual 
analysis comparison was also used in determining that the 
Pomme modes had been detected in either direction. No dis- 
Srepancies were noted in any of the systems checked on the 
masis of these comparisons. 

Program VIBREL includes in its output a graph of 
frequency determinant versus frequency in order to check 
Meetner the frequency increment iS SWech that : mode on 


modes have been skipped in the scanning process. =In tne 


om 








first analysis of a branched system which appears in 
Appendix F, the trend of the output plot showed eine. < 
natural frequencies in close proximity had been overlooked. 
When the program was rerun with the starting frequency incre- 
iWeme decreased by a factor of 2,7 thesemmves! requcmetccencmc 
eemected. 

The maximum number of natural frequencies which can 
Demtound using VIBREL is related to the Significant figure 
Capacity of the computer and the physical configuration 
feerne piping system. AS pointed out by Kime [Ref. 6], 
pemeene higher frequencies the columns of the frequency 
deecrmimant approach the point of being parallel; hence the 
itaiermeal walue of the determinant includes difterences of 
large numbers. Systems with few COMpPONents getmemakily 
em bit this parallelism at a lower frequency than the more 
@emplex ones: however, there is no way of determining before- 
hand where numerical difficulties will be encountered. When 
memeallelism is approached, a scattering of the values of 
meeauency determinantewill be noticed in the output plot: 

Ailthewenethe exponent Limit on@iene IBMesG0” cempucer 
feed Lor testing the program is +78, program VIBREL has been 
memeac tO Cease Computation for a particular problem and print 
a message when the exponent of frequency determinant exceeds 


+60. 
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D. CONCLUSIONS AND RECOMMENDATIONS 

Based on the reswits of the accuracy gangs te tna, 
analyses previously discussed and through comparison with 
Similar programs developed for planar cases, namely those 
@oerink [Ref 5] and Kim |[Ref. 6], programe PR Ric emlamne 
aiployed in its present form as a workenes tools 1omecnaimce n= 
in the dynamic analysis of spatial piping systems. The 
meer tiat Cighteen mode frequencies were obtained for the 
complex system exhibited in Appendix F implies that at 
Meast the first few and usually the most usefull frequencies 
mame determined for practical piping Configuration with 
mmeaecuracy consistent with engineering design. 

Pinor MOdtitcations son sadapringmthe program coed 
Pemticular situation can be made using the guidelines of 
Pyemendix C. With regard to major revisions, program VIBREL 
was developed with solution accuracy and integrity of first 
importance. It is recommended for future users that changes 
Mmimcne coding be made for a minimum computer execution time 
Poaiat becomes a prime requisite. Avenues of Investigation 
meearding numerical difficulties at higher frequencies are 
also open to the future user. A similar analysis for the 
planar case was undertaken by Kim [Ref. 6], who utilized 
alternatives to the frequency determinant approach. 

Pltchowen ene caleulatiome of MOdcueiancsc OL Ene piprle 
System has not been included in the program because of time 
miittations, an outline of how to ine@mperate this feature 


is discussed im the following paragraphs. This modification 
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to the program 1S.recommended (45 4) further siinme acai tee 
in the dynamic analysis of 3-D piping systems. 

Let us consider again, for simplicity, the two-dimensional 
system of Fig. 2:3 whieh 1S Yvepeated below ie toundetnat 


me total system trequency Condition could) be sexmnmece amar 


= 0 (Ba) 





anes 45 Se Salvealls geeeehoely i ieleelie Wees 1. Sin 


Peeemashing determinant of the coefficient matrix yields a 
natural frequency. At this frequency there is a non-trivial 
Berution to Eqs. 3-1 and A,/A, GalnepDe st Ould. None lt cen 

A, and A, gives the modal values of the non-zero quantities 


Gumene left end state vector, A,* and A,*. Si ee tS inlays 


iaeezing condition would be At + AS = 1. With these values 
known we proceed along the main member calculating the modal 
Memlections and forces at each divisional point. This is 
done simply by multiplying the normalized and compressed 

mec end state vector by the system transfer matrix (2r x r). 


in Fig. 2.3 the modal forces and deflections to the left of 


point 1 would be 
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jee jE! 
W b b A.* 
_ [P22 P24 2 eo 
Vee | see 
M b b 
Z 1b 43 44 


Ingorder to permit unequivocal 1nGenereta lion lor mp 1 colar 
mmepit, it would be advisable te incorporate a point number- 
mieescheme tor each point at which modal forces and detlec- 
mons are desired, In addition, these forces and deflections 
weil be expressed in the local coordinate systems and a method 
mmouid be incorporated for transforming them to the global 
System, 

Bien a branch point, suci as Vein Pipe 255, 1s cncolmecwcar 


mitesuranch transfér matrix up to point 1 is formed. Thus 


W c Cc B 


13 14 it 
W C e B 
_ 25 24 2 (523) 
vo See 
M C 
ZS1p as 44 


where By and B, Bets! 16 lS Mbp colonrl OMeelieices 1 le Sieel ee 
meetor at the remote end of the branch. Since we know that 


mieederlections are continuous across the junction point 


‘t : | (Bat) 


Sip ae 
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The quantities at 1L are given in terms of the known values 
A,* and A,* so that the values of By and B, which are intrin- 
sic to the 1B deflections may be solved for through the 
relations 3-4. The known values of By and B. may now be 
used to find modal forces and deflections in a manner Simi- 
lar to that used for the main member. 

yust to the rightsof. the june Gili Oe  eeecus cag 
deflections are the same as at the left; however, the forces 


fee be calculated by adding the IL forces to the 1B forces 


mimcne global coordinate system. 
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APPEND 


NOTATION AND NOMENCLATURE 


The following Matrix NOTA tee meee s ccleeie Olle misc te 


{ ] - matrix 
{} = nh xX eolumnm veeuen 
DIAG { } - square matrix having the given elements 


Iie. tS prameipal diagenalsandwezero 
elsewhere 
— SUPEYTSCHIPE CenG GaiCelio MeGanispocemot a Orvon 
matrix 
ane fOllowine VeeronenOtatloOnmails 9UScCGmm imine mt c x Ur 
A - vector A 


-_ 


[A] 


lenig timed Wee tom. 


>t 


-B - dot predictwon vectorsm and B 


-_ 


x B --CrOssmpl Cduct. OF VeCrOrs anaes 


1 


ible A, l first Wists She symbois msed sim crener stnc 
fee Or Appendix B, with the corresponding computer program 
tmemerable identifiers. Following these appear the program 


femercable names which have no counterpart in the text. 


TABLE 22% 
EXT PROGRAM 
SYMBOL VARIABLE DESCReE I TON 
A AREA CLOSsS SeCCtl ond agea ouapipe 
A A finst Of two Vecuonemuc. imo tas pupae 


dimectional chance 
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TABLE A.1 (Continued) 


TEXT PROGRAM 


SYMBOL VARIABLE DESCRIP 
B translation matrix 
B B second of two vectors defining a piping 


directional change 


C compliance matrix 

D VECtOr Of Cenenalizediedies hace ieimnus 

an UX Uiae Vector fn locale da, ecernrem 

IE Jepay Young's modulus 

F vector of generalized forces 

G G Shear modulus 

GC G2 Drdamch coordimate transformation matrix 

Shy G3 branch coordinate transformation matrix 

Gul Gl branch coordinate transformation matrix 

i Nee sieroieee le oe Zoe =olily Se Lem 

DIX mass moment of inertia of a curved pipe 
about X-axis 

ins inertia tensor about point © 

I, Gore dale nimaG y. aici i a 

1 Padus Or Cyd 1Ommos pace  6mOse sce milo 

j Die VecCrOne ime locals. dirccrilon 

J Shy) Mone htanOt ecm iano panier se Gildon 

k unit vector in local z-direction 

k BFF bend flexibility factor 

Ni DL lencth Of a2 Stratoiencecerion of pipe 

L POtaAcLOn Mmaerix 

m DM mass of a section of pipe 








ie XT 


SYMBOL 


— 


PROGRAM 
VARIABLE 


UX 
UXP 
UXDP 


U, Ul 


TABLE le CCOnmeimic d) 


DESCRIP ILO 
moment about local y-axis 
moment about local z-axis 
axial force in local x-direction 
Gm x6 NU 1 sina ees 
point matrix 
GaGiIUsS Of CMG a Pll Calor s nlp Cem ind 
mean radiusmor pipemsect ren 
Submatrix of branch transfer matrix 
submatrix of branch transfer matrix 


Tadius of ,center Of mass Of a curved 
SCIELO Oe TIES 


Nal Menbemy Una taV eC imag s Gig Ol ato 
Draneiy junetron 


pipe wall thickness 

torque about local x-axis 
deflection in x-direction 

unit dyadic 

Unt Vector am local x -cqivectron 
Unit» Vecror iInedlocal x ’=direction 
Ummeeviecceor Inmlocaien odincerlon 
Emeans fer siamese 

Gderlection in loecaiey direction 
field matrix 

Shear force in local y-direction 


Smear force in Poca -ailrecticn 
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PROGRAM 


SYMBOL VARIABLE 
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<-¢ s+ DO 


AL 
ALI 


GAM 


UU 
Ins 
ODF 


RHOM 


D2 


TABINE VA. 1 | (eont ined ) 


DESC UP TLON 


deflection in local z-direction 

Sieaiee VecLor 

x-axis coordinate transformation angle 
x-axis coordinate transformation angle 
Z-axis coordinate transformation angle 
frequency determinant 

angular deflection about z-axis 

jWekoviswe joleue Wdame leeneicley tame Moll jo\s 

Poisson's ratio 

Shear distribution factor 

mass density 

angular deflection about x-axis 
angular deflection about y-axis 

Cire rar shire amen y 

array of section properties 

arc length of curved pipe 


distance from bend point to tangent 
JOO) Lote 


he. x 12> xX. aoe OO Gematc ee rans Omilcac Tom 


Naser 1 x 


lee le Z2Z-aAxloecoorednmate transtorpmartionm 


matrix 
value of frequency determinant 
outside pipe diameter 


Mmside pipe dilamocer 


See 








ESOT PROGRAM 
SYMBOL VARIABLE 


DD 


FN 


NID 


NIT 


NIV 


NK 


NMO 


TABLE A.1 (Continued) 


DES CRT inert 
frequency determinant eamray 
OGiay Oe natin a lai Regul memes 
gravitational constant 
SSGtIoOn dent itTeny code 
branch boundary condition code 
boundary condition sequence number 
Dbeaneh repeat daseraineame 
left boundary condition code-main member 
Diane Onl Gadel i? tlet mcs em Uinli Tamit 


dimensStonees "auray Of scaleme, dete nian amr 
values 


right boundary condition code-main member 
curved subsection override discriminant 
sitear deflecetomy rotary inertia disenimameame 
Stralgnt Seeelon tes ue a1ser imi ame 

See b VOMe Mum 

WeECtTOr OF HOUNdany Condition codes 


Mt ea Linke NCE LO 1a lS Ci Inne | totem 
branch 


OC WVE Te eMC ei Cece 


MUMbDor TO etter lomomrcquared fOr Miodmer 3s 
method 


Tihelalounit VeeGbtonmmasS ec rama: ot Om. nee 
member 


dais Dranch discrimi aine 


number of modes requested in analysis 
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TABLE A! a Con tamned:) 


TEXT PROGRAM 


SYMBOL VARIABLE DESCRIPTION 


NN iteration S€cCt lon enimbpen 

NPROB number of systems to be analyzed by 
program 

NS ie section number 

NSS mumber Of CUrVed pipemsubcecc1oms 

PX X-coordinate of a working point 

PY Y-coordinate of a working point 

IPod Z-coordinate of a working point 

RHO weight density 

UA POdt Cana x eto metia eC nne 

UB UC Y Pedi tee Or ed eli c Cusp sape 

UST Straight pipe transfer matrix 

WINCR me Clone yeemeine a eiile Int 

WN Maer ale? reqteney, 

X array or frequencies fer OUtpUe oles 

x array of frequency determinant values 


for output plot 
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APPENDIX -B 
CATALOG OF TRANSFER MATRICES USED IN THE PROGRAM 
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NE REND Tes 


DESCRIP LON, OP SE ROGRAR 


1. General Remarks 

Program VIBREL is a double precision FORTRAN language 
digital computer program which is capable of performing a 
meequency analysis of a three-dimensional piping system, 

Meee program accomplishes this analysis by means of the 

meister Matmux method an which e€aeh Sseetion ofepipe 1smeharac- 
momized by a matrix known as a transfer matrix. These matrices 
Wipe are functions of frequency are multiplied together 
mmomceession to form a transfer matrix which 1s characteris- 
mueeot the entire system. To combine Satisfactory accuracy 
with minimum computer time, a lumped mass model is used for 
meved sections of pipe while for straight pipe, a distributed 
mass approach is employed. Shear deflection and rotary 

iilemiia aré Optional in the model. 

Hie SOlUtION 1S obtained by IMm€rementing =the ws redquene., 
mi@miorming the system transfer matrix and frequency determ) - 
muemwith €ach anmcrement. The sign of the frequency determ:- 
ame 1S used to control the Search for and convergence to 
weecm Natural frequencies. The convergence acceptability 
m@emeeerion may be set arbitrarily by the user depending 
Smeethe accuracy desired. (See section 3.b of this appendix.) 
Peaponce WOULd E€XpDeCt, the greater the accuracy, the more 


wemoucer Cime required. 


oe 








The system is defined by working points which are 
read into the program with the intensive sand cxtemcayve 
mmoperties of the Section Of pipe pmecedimns ea ce litem ele veel 
for this procedure are contained in Appendix E> [me anten— 
eve and extensive properties may Not Wary Alona esascime le 
section of pipe but may vary from section to section. 
baopPective complexities, including branches emanating from 
Branches, cannot be handled directly by the program: 

ieee DOMNdaryweomdd tlONSen (ccm cc mot ban land Oeniee 
ieee Specitied by coded input for the ends of the main 
iember and the remote ends of the branches. Alternative 


Em@a@econditions are discussed in section 5.b of this appendix, 


4 


meee rogram Structure 
a. Main Program 

Pers ipl Veda row )Cliage tyes Sion me@atiie 10 6. 1 nee 
Mreseussion will pertain to that diagram. 

Vata whieh eontains working points. intensi yesand 
Seeensive properties, and boundary condition codes is read 
mMimemene program aS 1t 1S required for Computation. Local 
@eemarnate systems are set up at the starting point, prior 
momoand following each directional change, and at branch 
memetion points. From the working point locations and unit 
weenors describing the local coordinate systems, section 
memechis and coordinate rotation angles are computed. 

For computational purposes piping bends are sub- 


meectioned based on the length-to-diameter ratio in the following 
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Geometr 







pee Common Properties 
otorage and 
Geometry 


Seu, starting 
Frequency and 
select Frequency 

Increment 

























Form 
trenster ond 
Coordinate 
Transformation 
Matrices 


Increment 
Frequency 
until Sign 
Changes 


Combine and 
Calculate 
Frequency 

Determinant 








iterate 16 


_ a Pina curve 


4 zero of 
Frequency 
Determinant 


Output 
and 
Prot 


Figure C.1: Simplified Flow Diagram for Main Program 


Se 








Manner. 


Onculiy el 1 SUBSECTION 
1 = 60/70 oes 2 SUBSECTIONS 
Me Ln << 3 SUBSECTIONS 


iio /D 1s greater than six, the number of subsections as 
computed as twice the number of mode frequencies sought, up 
to a maximum of twelve. The user may, however, choose to 
override the number of subsections as determined by the 
peeeram through the use of the parameter MCS described in 
Ppependix £. 

As the lengths, subsections, and rotation angles 
are computed, they are stored, together with intensive and 
Seaeensive properties, in an array. This array is passed 
mtcecommon storage (dotted line in Fig. €.1) to funetion 
Ser outine FRDET (FRequency DETerminant) which later forms 
and evaluates the frequency determinant for a given value 
Pie requency. 

When the geometry for the entire system has been 
Gomemited, frequency iteration begins. The starting frequency 
meoreoitrarily set at 0.1 radians per second. To compute 
the starting frequency increment, the program constructs 
aesynthetic straight pipe, with a length eaual to the 
total length of the main member, and having outside diameter, 
wall thickness, and intensive properties equal in magnitude 
to a weighted average of the main member sections. One 


wean Of the fundamental frequency Of this Synthetic pipe 
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with fixed ends is taken as the starting frequency 

increment. The system transfer matrix and frequency deter- 
minant are formed by function subroutine FRDET for each 

new frequency. When a sign change is detected in the 
meequency determinant, an iteration process, Using imellen s 
memmod Of Successive bisections and inverse, pagabolicsinter- 
Moelation is used for convergence to the natural frequency. 
Mi@emeconvergence acceptability criterion 1s taken as the start- 
mye frequency increment divided by ten thousand. 

After a root has been located, the search begins 
amew with the starting frequency equal to the value of the 
natural frequency plus one radian per second. When the second 
memiemirequency has been located, the trequency increment 15 
Changed to one tenth the difference between the previous 
mrommacural frequencies. 

INNS Seehieele I OeCAss Cemiewimles imei inves See ws: iecl 
number of frequencies has been found or until the significant 
figure capacity of the computer has been exceeded. Then the 
mode frequencies are printed and the pairs of frequencies 
and frequency determinant values used in the search are 
sorted mid Dice sed 

ie otbroutines 

The program weilizes sixteen subroutines in the 
Emieess Of computation. A brief deseription tollows: 

ANGLE Computes the sign and the angular 


Sepauacion Of Evomeemer Or munl te syvCcctons 
which have a common axis. 
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COORDX 


COORDZ 


CURMAT 


BELER 
DROOT 
ENV ERT 
MATMUL 


EOUNT 


ERPLOT 


SORT 
STATEM 


STRMAT 


HUVEC 


Forms an n X n (n is some multiple of 
5 and < 12) coordinate tranctormatnem 
Matrix for a rotation cabemesaelocar 
X-axis. 


FOrmMS an n xX nN COOTrdinere trans toma ton 
Matrix for a rotatione abel taeeloe alae 
axis. 


Constructs the fie ldgmatriose one weurnvecd 
SC Gul Olio t stn ee 


Computes the value of ann x n determi- 
Nant; Wsed to evaluate mene ot Peaienc, 
Ge te Vii antee 


Performs successive bisections and inverse 


HaGavolreyineempelatilone tOmlkeGaiccs thie 
ZETLOS FOr ete at Veg Ueney dete mmlinan tt ate a 
a sign change has been detected. 


INVERES. Cie unatr i> Ras Losec  loeadain 
SCompuring the Dranchepoumt una tiie 


Multiplies two matrices together. One 
of the matrices must be square. 


GCOMSELUGtES 4 soot matt 190) -Onesae laced 
SECtlon: OF epic 


Takes the array of frequencies and 
frequency determinant waluces veqlleu lated 
during iteration and plots them length- 
Walsis Psbiells eyo) lionesse Seely 2S Gio 
au cOnaticalily sand tthe Hemera sp let 
may be changed to suit the user, 


Serts the array of ireaquencies and 
frequency determinant values from 
ie Gal bal Ones Of ep Ot © imloe one le ssp Tac ln 


Constructs the state matrix for the 
Sp OrOD Eid te dase der tc cmmmOk =O lle Omit 
boundary conditions. 


Constructs the teancierematcrix Lor a2 
Straacht Secticmuomepine. 


Computes a set of orthogonal unit vectors 
from two pipe section vectors which 
Fepresemt a directionemeer mec in the 


jet Olea 
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WSTART Computes the value of the frequency 
determinant for a previously constructed 
Synthetic pipe equal in length te the main 
member and having weighted average values 
of the properties of the main member 
sections. This frequency value is used 
in selecting the stammimc 1 peauene, 
increment. 
ERUET For a particular value of frequency ,; 
this subroutine assembles the state 
matrices, transfer matrices for each 
SCeEtIONn Of PLDC we aiGeeeiemG© Chaentchiac 
transformation matrices. It then 
multiplies them together in sequence 
and calculates the frequency determinant. 
An interesting phenomenon occurs in the formation of 
micmiunrequency determinant for a single straight pipe. 
Meese Of the symmetry of a pipe in the y and z directions, 
the two flexural submatrices in the straight pipe transfer 
Matrix are identical, which causes the frequency determinant 
tO maintain the same sign on either side of the flexural 
Natural frequencies. This isolated case is handled in sub- 
routine STRMAT through the use of the straight section test 


aqiscriminant and appropriate exclusion of one of the flexural 


Suomacrices. 


Seeeremarks on Instructions for Program Use 
meee General Remarks 
Detailed instructions “for program use are Contained 
in the listing of Appendix E. It was intended for these 
imstructions to be part of the program in order that any 
Meer, given only the program listing or the computer card 
deck, could employ the VIBREL package with no prior knowledge 


eameranster Matrices and without perusalee: thisethesis. 
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Some claritficatwen of the instructions sone pree aman 
use 1S justified forMeases where properties  oneuceneim:, 
Seange within a pipime bend, First, @onsider spices ioe 


ceemmon of piping showm in Fig. C.2eqASsumiie That no change 


a 





mien. 2:  Papning —Renduwith No Giawee iieProncrties 


wae oss Section or properties occurs between points a and 
Geeeene only points which must be listed as working points 
imp rogram input are a, o, and d. Using these points and 
miemradius of curvature, program VIBREL will calculate 
Hemtieernis ab and cd as well as the arc length and anelle of the 
bends. 

Now.consider the scctaon Of pilpangeet erie erm aieme 


Micmeross Section or properties change at point c but remain 





e 
Puourewe 15: bipane Bend with wGisaqemo um monet eircs 
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constant between points a and c and between points c and e. 
Inethis instance points a, Os O95 and e must be listed as 
working points for the program. VIBREL calculates @stranene 
lengths ab and de as well as the information required for 
curved lengths bc and cd, It should be obvious from this 
example that for the general case of two curved sections occur- 
ring sequentially that the intersection point (point c in 
Fig. C.3) need not be listed. 

If the geometry input is such that negative Jengths 
are computed, an error message iS printed. One half inch 
Meemay 15 allowed for negative lengths to take into account 
any errors which may occur in measurement. 

eae Ycoeramn Fhexibilaties 
(1) Approximating Small Piping Accessories 
small valves, flanges, couplings and other piping 
accessories whose center of mass is relatively close to the 
centroidal axis of the pipe can be approximated using a 
penaehnt Section of pipe. The length, mass density, diameter 
aneewall thickness can be artificially varied to give the 
Maes ot the item. The resulting variation in stiffness 
Memveen an actual flange and artificial pipe, for example, 
tomeeeiy will not have a significant effect on the results 
Beeause of its short length. 
(2) Approximating Alternative End Conditions 
When the ideal boundary conditions are applied 
memene State vector on the end of a system, it must have 


eee eZero and six non-zero elements. Actual field conditions 


Os) 








may dictate that this 2s mot the Case because sore: leven 
mountings or connection to machinery. A good approximation 
in these circumstances may be made by attaching one end of 
am artificial Straight pipe to the system oneprmg scope ie 
omer end fixed, The mass of the artifvewal pipes nemidape 
small with the extensive properties and elastic modulus 
mileacn tO approximate the stiffness of the mounting. 

Thespregran Can then usemene stdea ls poumda iy, 
eemearcions with little loss of accuracy from the real situation: 

(Ss) ePreer aneiodi+ lcatroms 

In some instances, a user may desire to change 
Memeduantities as frequency increment or input format fo 
eweesnisS particular needs. For this reason, Table C.1 
lists the more important program variables and what card or 


cards to make these changes to. 


TAR ie oC 
ne : Card Variables 
MOE ELE RIEL Numbers Poet eae 

Maximum number of points for 2890 x(400), y(400) 
Pent plot 
iemimum number of piping sections ZOE AA(100,10) 
in analysis 
Maximum number of branches in 2920 NBCSG) 
analysis 
Input format 3210-50 
Output format (Input Data) 3390-420 
Output format (Properties and PES de od 
Geometry ) 
Oieeput format (Natural 8740 


Frequencies ) 
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TABLES CeigetComerlniued) 


Modit vcativon 


"eorttineg £requency  Inerement 
Bieanrting frequency 


Gomercence acceptability 
eraterion 


Maximum iterations to 
Gonvergence in DROOT 


Sieertinge Frequency after 
Mgeation of natural frequency 


mecduency increment after 
second mode 


ime r Of lines of print plot 
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Ganac 
Number 


8650 
8760 
OZ e 


3670 


9330 


9320 


14730 


Variables 
Affected 


WINCR 
W1 


BP S 


hee 


W1 


WINCR 


RY 








APPEND LAD 


PROGRAM FLOW DIAGRAMS 


ieee General Remarks 

Flow diagrams are included for the main program and 
mime taon subroutine PRDET. Subroutines DETER DROOT, 
MWeERT, and PRPLOT are standard package subroutines with 
minor modifications made to accommodate them to program 
feo. Since their equivalent or an alternative routine 
mayebe substututed for them, flow diagrams are not included, 
Adequate comment statements are, however, contained in these 
‘moroutines to define their structure. @ther subroutines 
fomewihich no flow diagrams appear, incorporate, at most, two 


@eeisions in their structure, - 
fee iain Program Flow Diagram 


Pages 6f through 76. 


Seeeeunetion Subroutine FRDEY Flow Daagram 


mages 7/7 through 79. 
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\ WRITE: 
HEA OING , 
INPUT 


NB (1) = Be 
NB (2)2 MES 





\REAO: Mea, 618, 
ere fee, oe 
Qty Qa, En 





—- --- 





' 80) * Px8-Pxi2) BC)= AG) 
8(2)*AZ) 
863)* AB)+1.0 







| B(2): P18 -P¥(2) 
803)" P28 -F2(2) 


ae 











AQ) PYC2)- PYG 


Ati) * Px(2)-Pxd). 
A(3)= Pe} Be ‘) 






CA, 8, UX, UY, V2Z,GA4} 


: 
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| oi oleae (7a ieee 

YRS cian Spat At : ia | 


<b Bp mm nh ee OED 
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Ne 
Se 
Oc: DSQRT (Ab) ve 
E A(2)K #2 4+ AZ NG, 


















AA oe 4)= 4.00 
AA(N2)* EY 
AA (1.3) = FE 
AA (N,4)= RHO 
AA(W.S5)* OL 
AA (NG)= T 
AAW,7)= DL 


END f0jNT 
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24 = OSQeT 
CAG) x2 +A 2} 
+ A(3) "#2 
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AAW,1)? 400 
AA(NW_2): EY 
AA (N,3)* FL 
AAI(N,4)= RAO 
AA(N, 5)? OL 
AA(N,6)= T 
AA(N,7)= OL 







Px(t)* PX (2) 
Py): PY(2) | 
P2(1)= PZ (2) 


QHANGE OF 
PROPERTIES FOINT 


OL=0O4-8L 





— 












es 
BEND N oe ue 
POINT >| 
BU) = PxB - Px(2) [| ux PCI) = Ux(T 
| BG2)- PYB-PY(2) | uyP(T): uy ae 








eae vz oS 
eer 


Chil UVES 
A, B, x P,UYPUZFPG AM, = 
ALPHA = 9.00 


vsctase a 
MOY, UY PO BPH | 


| 803): £28 - PZ) | 





ARCL = ket Gh 
641 = Be 













V2(L)= UZPCL) 
UX (TL) = UXPO) 4 OL GAM) = OVP OQ) * ASIN (GL) 
Uy (1) = 2X P(D)« OSIN GAM , + OVP(D) 4 LOS EFM) 


iy DELTA © F7- Gan 
|BL° R/OTAN (OLLTH/2 
KA*ARCL/OL 









ENE WEY ics 
Put OSAKT 
(AG) xz + AZ HER 
a ; 4 AC3)K*x 2) 
AA 1)7 2. BAN 6) = Ta 
AAN2)* EY RANI): A 
Bea) FE AA(N2)* Ark 
AA(N4) 442 AAW1)* GHA) 
AA(NS)* OL ABN IC)* OF WAT 
(NSS) 


a wes - —— 





AEE een 
N= Net | 





i) ae AA(N,1)* /. AAW, 4): LAO 
Px0)=Pxlz2)  VEY=QEY | PAG Ey AA(NS)= Ot aie 
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APPEND Tor 
PROGRAM ACCURACY AND INTEGRITY ANALYSIS 


1. General Remarks 

This appendix tabulates by mode frequency (radians/second) 
mnie results of the accuracy and Snbecnme, wanalyces) of -vaimmena 
piping systems from simple to complex. 

The comparison frequency values for the simple Straight 
Section systems were obtained by solution of the differential 
equations for longitudinal, torsional and flexural vibrations 
with application of the appropriate boundary conditions. These 
molutions were evaluated on the digital computer with accuracy 
meat least as many places as listed in the tabulation. Evalu- 
meron Of the comparison Solution and VIBREE solution were 
fried Out using a value of 32.174 nDeGe for the gravita- 
tional constant and a value of m7 to US sionviicantetieiigess 

Comparison values for the curved pipe analysis were obtained 
from three sources: Fink [5], Kim [6], and Ojalvo [8]. Those 
extracted directly from [5] and [6] were compiled Using a lumped 
mass model and digital computation. From graphs and formulas 
(based on the theory of clamped ring segments) appearing in [8], 
mintemouhner natural frequency values for the curved pipe were 
Galculated. 

Beveral comparison frequencies for agsinple braneh System 
were obtained from [6]. 

shear deflection and rotary Iiertia were neclected jexcepe 


where noted. 
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le 


SYSTEMS ANALYSIS 


a oFRAIGHT PIPE With VAR LOUS 


—————} 


END CONDITIONS 


SYSTEM PARAMETERS: 


DIAMETER (IN) AW 
. WALL THICKNESS (IN) MIM 
END CONDITIONS: LENGTH (IN) ays 
FIXED- FIXED MODULUS OF ELASTICITY (PST) 
POISSON'S RATIO : 25 
WEIGHT DENSITY (LB/FT”) 490. 
COMPARISON RESULTS: 
MODE ANALYTICAL VIBREL * DIFFERENCE 
fot FLEXURAL 75.10461420 75.10461590 .00000 
2ND FLEXURAL 207.02876110 207.02876110 .00000 
3RD FLEXURAL 405.85914940 405.85914210 .00000 
4TH FLEXURAL 670.90578980 670.90577600 .00000 
5TH FLEXUPAL 1002.21750380  1002.21750410 .00000 
6TH FLEXURAL 1399.79137910 1399.79140510 .0n000 
7TH FLEXURAL 1863.62757630 1863.62799710 .00002 
1ST TORSIONAL 2007.8325054 2007.8275305 .00025 
8TH FLEXURAL 2393.7260868 2393.7157896 00043 
OTH FLEXURAL 2990.0869112 2990. 3653337 .00931 
1ST LONGITUDINAL 3174.6619386 3174.6560265 .00019 
DUAL ANALYSIS RFSULTS: 
STARTING AT STARTING AT 
MODE LEFT END RIGHT END * DIFFERENCE 
1 75.1046142 75.1046142 .00000 
2 207.0287611 207.0287611 .00000 
3 405.8591494 405.8591494 .00000 
4 670.9057898 670.9057898 .00000 
5 1002.2175038 1002.2175038 .00000 
6 1399,7913791 1399.7913791 .00000 
7 a 1863.6275763 1863.6275763 .00000 


as. 





4 ________— 


END CONDITIONS: 
FeexeED- FREE 


SeMPARISON RESULTS: 
MODE ANALYTICAL 


ci 2 6.657 
OG STi 
Zl Gare 
F058 5409057 
6) JOE SHO eRU NG iS) J 
eo 7 OL SFOS 


Am WYO 


mUAL ANALYSIS RESULTS: 


STARTING AT 
MODE Ter TEN 


oz OO 
(ENG NS (oS 
Zell OG 
aI )S) 7.10) Sel SHS, 
OU 0 oUie Sl 
INS) SHS si oyestee 


NM PWN 


SYSTEM PARAMETERS: 
SAME AS (1) 


MO ievecte lb 


LILA SOZ S606 
(POP so) SVS) 
PUY BRANES SS) 7) 
405765455 
6070. 99G0 G51 
139 Oe desis s 


STARTING AT 
ReG rite fei) 


bre Rees, 
1560107 S057 
ZO] Fela Orc, 
403. 8541956 
670 900 Gios1 
N39 9.7 SMES 5.5 


IAS 


DDPPERE NG 


.00003 
.00002 
mIOO 016 
.00000 
SOO 0 
.00000 


DOE PE Re eis 


.00000 
.00000 
.00000 
.00000 
.00000 
.00000 





<>, 


SYSTEM PARAMETERS: 


SAME AS (1) 


En CONDITIONS: 
fet. D- PINNED 


BeMPARISON RESULTS: 
MODE ANALYTICAL 


Seo) LOUD 
167.7264474 
549 .9478448 
Bytes 4 1S) ay 7at0 | 
ios | he 
MOUS 0 Log S70) 


Nm PWN 


MUA ANALYSIS RESULTS: 


STARTING AT 
MODE JE It ed, 


alee? s/f 
167.7264180 
SA a4 78400 
DO od O22 6 
Sr IN arses iy 0 
OO JE Oey 0 


ON UT SND Fe 


VIBREL 


Dele) Ore wes 
167.7264180 
349 .9478460 
oor, 45 2 Ao 
O15 24 7 7 5 oan 
OOS 2 NO io 


STARTING AT 
RIGHT END 


One. LS ieee 
Voy. 7 Zoe w 
549 9475400 
598 44 Soa 
OA. 7 one 
1003. 9164796 


ee) 


DIPPER ENCE 


. 00003 
.00002 
00000 
~ 00000 
.00000 
SUC Ue? 


DERE er 


.00000 
. 00000 
OOOO 
. 00000 
OOOO) 
.00000 





N 
N 
(-35,3 - 35.3, o) at 


— ae ee ee eee eee ee ee 






END CONDITIONS: 


SrolEM PARAMETERS: 


~ (120.6,50,0) 
/ 


(35.2, - 35.3,0) 


JID 12 Eu. 


PIPE DIAMETER (IN) ae 
WALL THICKNESS (IN) ive 
INCLUDED ANGLE OF ARC (DEG) 270 
MODULUS OF ELASTICITY (PSI) 30 x 10° 
POISSON'S RATIO : LOS 
WEIGHT DENSITY (LB/FT”) 556 
COMPARISON RESULTS 
a. MODE 1 
SOURCE FREQUENCY 
VIBREL 37.2234058 
OJALVO [8] 38.5930983 
FINK [5] 38.8773520 
KIM [6] 35.3639740 
b. MODE 2 
SOURCE FREQUENCY 
VIBREL 72.1914615 
OJALVO 69.9962862 
FINK 73.3297307 
KIM 66.6993490 
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* VARIATION WITH 
VIBREL VALUE 


Ole 
4.44 
= Send) 


VARIATION WITH 
VIBREL VALUE 


one 
Ly aye 
EW hs OnE 





C. BRANCHED SYSTEMS 


(1) SINGLE BRANGH=S voc 


(29.3, 70.7,0) 





SYSTEM PAR AMET TRS. (Alias =o) © PEON Gum 





DIAMETER (IN) Pate 
4 WALL THICKNESS (IN) Phos 
(0,0,0) (100,0, oO) (200,0 6) LENGTH (IN) 100 
MODULUS OF ELASTICITY (PSI) 30 x 109 


POTSSON 'S RATIO 25 


WEIGHT DENSITY (LB/FT°) 556 
COMPARISON RESULTS: 
MODE KIM [6] VIBREL * DIFFERENCE 
2 194.25470558 194.2888499 .01758 
5 278.61078241 278.6390155 .01013 
8 628.59556355 628.9417782 .05508 


DUAL ANALYSIS RESULTS: 


STARTING AT STARTING AT 
MODE Peet END RIGHT END o URE PER ENG 
i Of 04 002 / 67 EO 2S Coc .00000 
Z 194.2888474 194 ..2886474 . 00000 
5 oe Oo es 55 LOS 730 Cae os 00000 
4 2A | a G9 242 0G .00000 
5 Ze SOO UES So 27 3-0 SSIES 3 . 00000 
6 281.9864461 28159802461 . 00000 
7 S70 (9492s 53/0... (9472s .00000 
8 OS, I fee Z O26. 2a aia .00000 
9 634.7478785 O34. (47a og . 00000 
10 OG Seo olen Oa0 5/5 72 mine . 00000 


Ae 





(2) DUAL BRANCH SYSTEM 


(40,30, 0) END CONDITIONS: 
MAIN MEMBER - FIXED-FIXED 
BRANCH #1  - PINNED 
BRANCH #2  - FREE 





(90,-20,0) 


Seo1TEM PARAMETERS: 


SECTION SOTTO SECTICH SECTION 


PARAMETER 1 2 3 4 
DIAMETER (IN) (0 20 2.0 70 
WALL THICKNESS (IN) £125 £125 12S 125 
LENGTH (IN) 050 ees A ee LL 
MOD. OF ELAST. (PSI) 30 x 10 30 x 10 30 x 10 30 x 10 
POISSON'S RATIO £30 30 30 ae 

WT. DENSITY (LB/FT93) 490.0 490.0 490.0 490.0 


Weed ANALYSIS RESULTS*: 


STARTING AT STARTING AT 
MODE LEE END RIGHT END > DIFFERENCE 
1 HenOtes 0 534] TAO. SOA 554 7 .00000 
Z 276.6928874 276.6928874 -00000 
5 Sk Ago OS S0 SAS Soe 6 ~ 00000 
4 MIU Nis SEG Sis. 417.0890 52 .00000 
5 AVS) eh ats) abe (5%) £63 <5 20 celeZ .00000 
6 675. 0605260869 6/550 occa .00000 


PeeoeeAR DEFLECTION AND ROTARY INERTIA ARE CONSIDERED IN THis 
ANALYSIS. 


+ COORDINATES SHOWN IN DIAGRAM ARE IN INCHES, 
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COMPLEX, Sys Ean 
END CONDITIONS: 


MAIN MEMBER - 


BRANCH 


(0,0,0) (20,9,0) 


.* 
o 


(35,0,0) 


(60, 


FIXES XED 


ae G58), 
(60 25,-25) 


26,40) 


meres ANALYSIS RESULTS: 


* MODE 


MOAOnNIAMPWNFE 


STARTING AT 
Pelt SEN 


Tone 
IO ar 
142. 
ZO. 
29 
yaw 
Sole 
364. 
Eto 
a2 
542. 
Seo 6 
eee 
648. 
Laie 


350449 

7404907 
5065200 
1255 ho9 


.1786608 
.§210841 


Selig 20 
ao Oa 
5022504 


~1734071 


AO GS OF 
0762600 
0701684 
8240390 
§640271 


~ 


MOS2 04560787 
IANS PAO ass 2) 
ZO E28 Go cn, 


(60,25,/0) 


R23.0 


(99,-9, 55) 


SYSTEM PARAMETERS 
(ALL SECTIONS): 


DIAMETER (IN) 2s 
WALL THICKNESS (IN)  .035 
ELASTIC MODULUS (PSI) 30 x 10° 


POISSON'S RAI TO 3 ~ 30 
WE. DENSITY i0EB/ PTs). 340 ees 
Cbd. 
(/4-0,20,-/0) 
222.0 
(//0, 20,20) a 
(/0,20, 20) 


(/00,-20,60) 


START Gre 


RIGHT END % DIFFERENCE 
73. / ota a0e OSS 
102. 9 2G so a Genl 
W46e SoU Ges 60 2.674 
201.3564078 es 
239.1641100 006 
62121 Geass PAPAYA 
O57. so cer a0 JL 80: 72 
SoA wis 5 oe IEG, 
429.1841008 145 
B22 one Ge 207 
54452 00m O25 
Bigg Ore 23 
57 0 70S 522 
6514 28oice2 CHES, 
759 2450 eo 606 
WOSZ Sez 0 026 
2 Oy Sze 016 
124 7 S520 cs 224 
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APPENDIX G 


SAMPLE PROBLEM 


ime Piping System 

The piping used in -thuseexanplesioedcmenead commen. 
Semplex system shown in Appendix F. It includes several Corners 
and bends and has one branch junction point. The boundary 
conditions are fixed for the two ends of the main member and 
Gee remote branch end. 
2. Data Cards 

Shown in Fiao, Gel US tEMe” anrancement of cauame tras 
as they would appear following the program deck. Note that 
Ward #l specifies one system to be analyzed, and Card #2 
metes) tOr four mode frequencies in the analysis. Since the 
Gemeyed SubSection override is not utilized, the program will 
go ahead and determine the number of curved pipe subdivisions 
meeerding to the length-to-diameter ratio, Shear deflection 


ama rotary inertia are not included in this problem analysis. 


GS 2. 













m iti, i ae: oh Oe el ee i es a mou x 
pn of =S.n ti me i not th. i een r,t a : 
a... - ; oe ‘ 
ihe SH. Snr oo nt oe Oe Oat is \ 
= Cl osl Ay i Ca | Fin “TT cis i ” am a ae ie 

Mm 4,4, ee fee Sig eee ee ie letra ee eet 

7, wo s oat a at ra a ee ae iaaaag ae: 

Meise Su il i BRANCH STA Ufo POPOV, 

da 2 a = e ag a i ee ot i 7 : 

at Pedic | — sol aya ch bt odin ft if Moi Hew 3 Tae eee . 

4 eaene: *. Sh pa! = a oe =) °, a =a >. = 

ae oe ae L en ee Roa ree 3. LTS 8 ee Si oa hisses 

ES a a a LS ee ny » a 

Ree a oy ieee i 
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i . ory hoyt 


{ - Se ee oF 
Ea 2 2777 a a 
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ri GAL 


mee, . 
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ore vt PR Weel hy a 
“oe a ; ‘“ F ‘Ve 
a we * ce ty ~ one LY as 
& ; ¢ ¢ ra - 
v7 LU '’ 


,00CTDCTHOVDCDIACDEGONDOODS MOTCC SPELL sere er TO) Ug 0UO NONI) Ce 


Sec eee aro “9 8 89 44 95 1G tT ae oO A 88 8 94 8 Or > oe te Sb fe A Oe Bt 5 Oo eR oy SR FE 620 See Lio Go (yr Sf coil tbiecua) a) viene soualec 


i isiiisitiicvcsecgeecasscrgugeditupeae Visage gat bi Melee 


M222 2222222222227 222222702: 


. 
Ld 
? 
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+2 
eo 
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rs 


1 a fi? 2792 2 22 22 2) ee? 2 ee 
emer ge 9 99999339333393333993I9532 


t.3 
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tw 
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tw? 
[oe 
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Cu 
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[a 
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MP SSAA LAAS ASA SAA SSSI CAAAE She LSE SS AMAA Sha cae esa daa 
§$55555555555555955 SSSSSSSSSET OT SRSSS USES SISESGSSSTSSSSESSSSSSSSSSESSSSSSSS ESSEC SS 
cs cecissiciale, (Recsvs elise viseesst SESEHSELESEEGSELESESSS 
PUTTTTTTVV2VTTTTTVVIVTVTVTT TTT LG TELE T TTT TTT TTT TTT 


— 


595992 393993 S38 so ee 


ee ciedieiiissaee 


29 : Gee 2 AO cre te 2 FG q wo 6 e0908 4 Ane aR GAs Que Guhed go ao 
Pires sassiccasesasessgasacasagaaaegsegsassz sag aess esa glee te} 9995559999899292 
ss ia ae Savi ea he 2c Serato Sea, Fiera tomate Sate 5d) 67049 34a) sede a cue ete es ot oN Gad Ne t et eet Seg: 


1EN 2345505 


= ss en a a ee ee EEE 


Figure G.1: Data Deck for Sample Program 
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RIGHT BC 
COOE 


CLEFT BC 
COOE 


* 
* 
* 

WALL 

THICK WT/VOL 


0.0. 


C.0.RUOOLF 
SRATION 
CLUOEO 


HRSA KKK AREER EKER EKER ERE EER ERE KEE EEE EKER EHH EH HERE H KH HH 


fa 


PROGRAM VIBREL 
COORO COCRO 


uf 


x 
COORO 


THE EFFECTS OF SHEAR DEFLECTION ANO ROTARY INERTIA ARE NOT CONSIOEREO IN THIS PROBLEM 
10 


KHKKEKHKH AERA SE KAKA KARE EERE EERE EKER AE EHEE EKA AERA E SR ERE HS EEE HK 
KEEEEAE SKK AREER ES SKEET EKER AEA REAR EAR EAE SSHKAAE ERA HES ERE ERE SES 


eee CCC CSTeTTeCTeT eC Ce TCCCCCSCTC SCO SS COSC CS ST CTT SST ee SPC CS eT eT ee Te Ce eee eee es S| 
Tee SEC SSS Se SESS STC SSCS IOS TSS OCT CCS SCC LCS CS STC CS TSS ISLES eS TCS SS ST Se ST ee eS 2 tS 


PROBLEM # 1 
INPUT OATAS 
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F CURV 


‘Ballo 
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WT/VOL 


PROPERTIES ANO GEOMETRY: 
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oO Se 
CoO0000000000 
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FREQUENCY OETERMINANT 


KKK AE HARE KE KAA KKK EKA KEARSE KARR HERE AEE EKA AEA ERE REHEARSE EKER KEE 
FREQUENCY(RAO/SEC) 
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PIPING SYSTEM NATURAL FREQUENCIES: 
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MOOE 


GRAPH GF THE VALUE OF FREQUENCY OETERMINANT VS. FREQUENCY 


1.270€ 


2.539€ 


3.809E 


5.079E 


6.349E 


7.616E 


8.886E 


1.OL€E 


1.143€ 


1.270€ 


1.397€ 


1.524E 


1.651€ 


1.778€ 


1.905E 


2.032E 


-5. 


ol 


Oi 


ol 


Ol 


Ol 


Ol 


01 


02 


02 


02 


02 


02 


02 


02 


02 


02 
-5. 


Wes -2.538E-17 3.585€-19 2-610E-17 5 »L84E-17? 
SHETEH SHS HELHEHS SO ESOH ETEK OHEE SOS OHH HHE KH OEKEHHHHHTS SESE SEKSHH RESTS EEK SKE LESSEE KEE SEE 
+ 
M so 
6 e 
6 6 
* e 
* 6 
6 * 
6 s 
6 s 
* e 
+ + 
6 A s 
* e 
s * 
6 e 
s e 
td e 
s * 
* ° 
td s 
+ + 
* * 
s 5 e 
s s 
s s 
* s 
* s 
* * 
s * 
s * 
+ + 
* s 
e a ° 
6 4 s 
s * 
s * 
* 6 
s * 
6 * 
* s 
¢ + 
* s 
6 * 
* * 
s € ® 
* * 
* * 
s * 
Ld * 
* s 
+ 2 
s r 
s s 
* * 
s * 
s 5 e 
s | * 
. L 2 
s A * 
* s 
+ * 
* e 
* s 
s e 
* s 
* * 
* i ¢ 
Ld * 
Ld * 
s A e 
s + 
* * 
* * 
s * 
* * 
Ld * 
Ld * 
* * 
s . 
* * 
+ - * 
s * 
e s 
* s 
e A ry 
* * 
cd] 5 e 
Ad Ci 
* . 
s se ry 
+ + 
s x Ls e 
* * 
s . 4 
Ld e 
Ld * 
s é e 
Ld * 
s a s 
Ld * 
+ 5 + 
LL * 
s s 
s A e 
s e 
s 5 e 
bd * 
s a s 
s * 
Ld * 
+ 4 + 
Ad * 
s * 
bd s 
s 5 ry 
s s 
s A s 
* * 
Ld * 
td é s 
+ + 
Ld s 
s e e 
* a 
s s 
* s 
s o s 
bd e 
Ie e 
* . s 
¥. + 
bd s 
* s 
s s 
* * 
s . s 
* s 
bd * 
s : ry 
Le e 
¢ + 
s * s 
* e 
® . 
@. s 
Ad e 
bd ® 
.., * 
Ld s 
* * 
+ + 
s Ps * 
Ld * 
* * 
bd . 6 
bd * 
iJ . 
6 a e 
* * 
« * 
+ . + 
SEE SEEPS STH HO EES EE EEH SHE SCHHK EES SO SEH EHEH OSE SESOHEREKRESS HEHEHE SHH SOHO ER EHEH HOEK SEE 
LI3E-17 -2.538€-17 3.585€-19 2.610E-17 5.184E-17 


LSS 


1.270E 


2.539€ 


3.809E 


5.O79E 


6.349E 


7.618E 


6.666E 


1.016€ 


1.143€ 


1.270€ 


1.397E 


1.524E 


1.-651€ 


1.7 76E 


1.905E 


2.2032E 


Ol 


01 


01 


ol 


01 


Ol 


ol 


02 


02 


02 


02 


02 


02 


02 


02 


02 


UU 


BIBILOGRAPHY 


Bishop, R. Es, 3. andy Johnson tite oc ase omen 
Vibration, Cmabridec wl naasemcs ty Pre > smc E 


Brock, J. Ese Complianeecmonr os) 1pe Components, | aNaggcal 


Engineers JoOUriialy. ses eeeoo. 7 WieOs ele od. 


Brock, J; ERgpediatrix Analisis cfertexabple FPalamemuse: 
Puoceedinegs Gumuie Hirst alo. Natomas Coneress som 


Applied Mechanics, p. 285-289, June 1951. 


Wen Hartoc, @ewra, “Mecchanteatl Vibpaetons, 46h Cds eects. 
fog, 1956% 


Fank, G. Ex 2vabration Anadlivsas of Piping Systems; 
M.S. Thesis, U. S. Naval Postgraduate School, Monterey, 
Calatornia, s964 . 


Loewy 2 One Virdee lOn Analisis Of PIpine Sys Cems. Mao. 
Thesis, U. S. Naval Postgraduate School, Monterey, 
Cali¢ormia, 1960, 


Vetracken, J, Do vand W. .- Dorn, Numerical Methods and 


Fortran Programming, John Wiley and Sons, 1964. 


Ojalvo, I. U. and M. Newman, “'Natural Frequencies of 
fmianped Ring Segments,” Machine Design, v2. SO. p21 o- 
eee, 2) May 1964. 


mestel, &, C, and FP. A. Leckie, Matrax Methods am 
Elastomechanics, McGraw-Hill, 1963. 


Bestel, &, Cy and F. A. beckie, “Transfter-Matrix Funda 


mentals,’’ International Journal of Mechanical Sciences, 
twee, Dp. 1o/sto/e LIU. 


L5G 





INITIAL DISTRIBUTION Liss 


No. 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 


Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, California 93940 


Professor John E. Brock, Code 59 Be 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93940 


i. . ebekeilea, ii 


2060 Sunrise Circle 
Marina, California 93933 


Lg 


Copies 


2 


INC LASo LE Tiny 


Secunty Classification 


7s GGUMEMINCONTROL DATAMMG IEE 


Security Classification of title, body ol abstract and indexing annotation must be entered when the overall report is classilied) 


ORIGINATING ACTIVITY (Corporate author) 2a. REPORT SECURITY CLASSIFICATION 
Maleiiasc ied 


Naval Postgraduate School aan 
Monterey, California 


REPORT TITLE 










Vibration Analysis of 3-D Piping via Transter Wetrtces 






OESCRIPTIVE NOTES (Type of report and inclusive dates) 
Master's Thesis; December 1971 


S AUTHOR(S) (First name, middle initial, lastname) 


Chester Davis Rudolf, IIl 


December 1971 50 0 


fa. CONTRACT OR GRANT NO. 9a. ORIGINATOR'S REPORT NUMBER(S) 





b. PROJECT NO 


c. 9b. OTHER REPORT NO(S) (Any other numbere that may bo assigned 
this report) 


10. OISTRIBUTION STATEMENT 


Approved Hom Ul eee ee ease ; distribution unlimited 


- SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 


Naval Postgraduate School 


Monterey, California 93940 





P43. ABSTRACT 


This thesis develops the theory for and presents a digital computer 
program capable of determining the natural frequencies of a three- 
'dimensional piping system having arbitrary configuration. The 
Migglysis uses the method of transfer matrices. Piping hangers, loops, 
fand complex branches (branches emanating from branches) have not | 
Mheen included in the analysis. A distributed mass model 1s Used 10% 
Meeraisit pipe sections, while mass is lumped for curved sections. 
1Inclusion of shear deflection and rotary inertia is optional. 


Several piping configurations are analyzed using the program, the 
fresults are compared with analytical solutions or values from the 
literature to demonstrate the accuracy and integrity of the program. 


= e 
CA A PS See rr = ate 2 as = - > - - ™ = 
at a ee ee le TTS vs t< age Ss a eNO oe Tags sa Se % ae 7 OO ity OS Tea 


DD 1 NOV wt4/73 NC LAS OLE Ree ee 


S/N 0101-807-6811 AR Secunty Classification A-31408 








N | al ee 
Security Classification 


eee 








meanster Matrix 
Piping 


Velprations 


: Pek Sia as - i “sb Sere ae! = 


We tas (BACK) 1 eg ee 


‘101-807-6871 139 Security Classification A-31409 














thesR835 
Vibration analysis of 3-D piping via tra 


DUDLEY KNOX LIBRARY 





